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Abstract 

O(a) QED radiative corrections to polarized neutral current leptoproduction of heavy 
quark flavours are calculated in the leading log approximation for HERMES, COMPASS, 
and HERA kinematics. Besides the inclusive case, we derive general ^-differential for- 
mulae for the polarized boson gluon fusion process and use them to calculate radiative 
corrections to semi-inclusive polarized heavy flavour production in the case of HERA 
kinematics. 



1 Introduction 



So far, the polarized gluon density Ag(x,Q 2 ) is nearly unknown, see e.g. In order to 
improve this status several experiments are designed to measure polarized heavy flavour 
production which is directly linked to Ag(x,Q 2 ) already in leading order (LO) QCD via 
the boson (photon) gluon fusion mechanism. Such measurements will be performed by, 
e.g., the HERMES Collaboration with its charm detector upgrade @ and the COMPASS 
Collaboration at CERN ||. Using fixed targets these efforts are restricted to rather 
low energies. A polarized HERA collider mode 0] would allow for measurements in a 
kinematical region not accessible to any other polarized experiment ||. 

Radiative corrections to polarized deep inelastic scattering using HERA and HER- 
MES kinematics have been recently calculated in complete 0(a) and in 0(a) leading 
log approximation (LLA) || disregarding heavy flavour production. However, at small 
x < 10 -3 the heavy (charm) quark contribution g^(x,Q 2 ) to the spin structure function 
g[(x, Q 2 ) becomes relevant. Keeping in mind that heavy flavour production processes are 
semi-inclusive it is recommendable to perform the radiative corrections on a more differ- 
ential level. As in the unpolarized case |7j we study the momentum (z D — p • pr>/p • q) 
distributions of the D-mesons produced in such events. For this purpose we derive the 
^-differential expressions (zh = p ■ Ph/p • <?) with general couplings, where h denotes a 
heavy (charm) quark, for the polarized boson gluon fusion process. This allows us to 
study electroweak (ew) radiative effects in heavy flavour production. Such effects have 
been recently investigated in the case of the light flavours u, d, s [|, |J and there turned 
out to be sizeable at small Bjorken-x (x = 10~ 3 , 1CT 2 ) p. 

Thus, in this paper we calculate radiative corrections to inclusive and semi-inclusive 
polarized heavy flavour production in LLA for HERMES, COMPASS, and HERA kine- 
matics. Since in leading order QED the polarized and unpolarized splitting functions are 
the same, i.e. AP ee (x) = P ee (x), the whole formalism can be adopted from the unpolar- 
ized case pj merely replacing the unpolarized Born cross sections by the corresponding 
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polarized expressions. 

The paper is organized as follows: In Sect. 2 the Born- level expressions for polarized 
deep inelastic scattering are given and the general ^-differential expressions for polarized 
heavy flavour production via the booson gluon fusion mechanism are derived. In Sect. 3 
we provide the theoretical framework to calculate inclusive and semi-inclusive (9(a)-QED 
corrections in LLA, followed by numerical results for the fully massive radiative corrections 
to heavy flavour production in Sect. 4. Finally, we summarize our results in Sect. 5. 



2 Born cross section 



In our calculations we will restrict on the longitudinal spin asymmetry dAa = (da(X, Sl) — 
da(X, — Sl))/2, where <ia(A, S^) denotes the differential Born cross section for deep in- 
elastic scattering of longitudinally polarized leptons off longitudinally polarized nucleons. 
For defmiteness, we consider the process l(k 1: \)P(p, S L ) — > l'(k 2 )X(p x ) where A and ki 
are the helicity and the 4-momentum of the incoming lepton, Sl and p the spin and the 
4-momentum of the nucleon and k 2 the 4-momentum of the final lepton. The Born cross 
section can be written as a contraction of a leptonic tensor L^ v with a hadronic tensor 

? C (Arr~i Z \ _L_ Arr ZZ 
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y xbXb> (iffWgJ) dxdy, (1) 



with B, B' — 7, Z or W , a the QED coupling constant, and x, y, Q 2 the well-known 
variables in deep inelastic scattering (DIS). xb collects some constants and the ratio of 
.B-propagator and photon-propagator: 

x,(Q 2 ) = i 

2 1 Q 2 

Xz{ - Q ' = (2sin0 w cos9 w ) 2 Q 2 + M 2 z (2) 
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where 9 W is the Weinberg angle and M z , M w are the masses of the Z° and W boson, 
respectively. 

Using a general 7 M (V — A7 5 )-structure for the leptonic current, L BB ' reads: 
L B J(\) = ^«(A,A; 1 )rf M (A / ,A; 2 )M(A / ,A; 2 )r^(A,A; 1 ) 

A' 

= 2Cf B ' (k lfl k 2u + fa v k 2iX - (fa ■ k 2 )g ta/ + 2i\e llvpa k{k 2 ) . (3) 

Tf = 7 M (V^ B - Af 75) and Cf B ' = (V t B - XpiAf)^ 3 ' - XpiAf) are built from the vector 
(V) and axialvector (A) couplings of the fermion-gauge boson vertex 1^(V B — A B ^ 5 ). 
They are given by 
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where Q f (Q e - = -l,Q e+ = +l,Q u , c , t = +2/3,Q d , s , b = -1/3) and Ij (I 3 e+ ^ t = 1/2, 
Jg_ dsb = —1/2) are the charge and the isospin of the fermion. pi accounts for the fact 
that the result for e + p-scattering can be obtained from the e~p-results by the replacement 
Af — > —Af, i.e. p e + = — 1 andp e - = 1. Furthermore, it is convenient to define two further 
combinations of the couplings at this place: 

Sf» = V f B V f B '±A f B A f B ', 

R B X = (V f B A f B '±V f B 'A f B )/2 (4) 
such that Ci = S t:+ — pi\2R t:+ . 

Assuming current conservation, the polarized part W M „(p, q, S) of the hadronic tensor 
depends in general on five structure functions gi(x,Q 2 ). In the literature different defi- 
nitions for the structure functions are used. For a comparison of several conventions see 



H]. We use the following form: 
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For longitudinally polarized nucleons the nucleon spin vector can be chosen in the nucleon 

rest frame to be S L = (0,0, 0,M). Neglecting terms of the order 0(M 2 /S) (S=CMS- 

energy) one can replace Sl by —XnP where A^ is the degree of the nucleon polarization 

PI and discard gi and g$. This yields 
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in agreement with the hadronic Tensor given in [11]. Using eq. (|I]), (§) and (§) one finds 
the following result for the longitudinal spin asymmetry (polarized part of the Born cross 
section) : 
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with the structure functions 

G 1 (x,Q 2 ) = \ J2 C ? B 'xb(Q 2 )xb'(Q 2 )9? B \x,Q 2 ) 

B,B' 

G 3 , 4 (x,Q 2 ) = Y, C i B, ^{Q 2 )XB>{Q 2 )9if{x,Q 2 ) . (8) 

B,B> 

The structure functions gf B> can be calculated within the quark parton model. The 
two relevant partonic processes are in leading order (LO) the polarized heavy flavour 
production via the boson-gluon fusion subprocess B*(q)g(k) — > Qi(pi)Q2(P2) and the 
electron-quark scattering B*q — > q. 



2.1 Polarized heavy flavour production 



The structure functions g^ =1 3 4 for a pair of heavy quarks Q1.2 with masses m 12 can be 
obtained by a convolution of the polarized gluon density Ag with the polarized structure 
functions g\ =1 3 4 of the partonic subprocess: 

" 1 dz? 



h.BB' 1 rfl 2 



h,BB' 1 r^2 2 
#4 " 



/02 2\ /" ^2 a / 2\ ~h,BB' 1 i rf2 2\ 

Q ) = / — Ag(z 2 ,n)g 13 [x/z2,Q ,n) 

J ax Z 2 

r\2 2\ f , 2\ '-h.BB'/ / ^2 2\ 

,Q,V)= — z 2 Ag(z 2 ,ii ) 4 ' {x/z 2 ,Q ,n) 

J ax Z 2 



(9) 



with a = 1 + (mi + m 2 ) 2 /Q 2 . As usual we identify the scales = /ip =: /i in Ag>(z 2 , A*f) 
and a s (fj,R) (see below). 



Following the general framework given in ref. |TT| , we calculate the polarized partonic 
structure functions gfl^f 4 in ^-differential form (h = c, b, t) by writing the two particle 
phase space as dR 2 = dz h /(8ii) with Zh = k ■ P2/k ■ q. The general z/j-dependence of the 



polarized structure functions turns out to be the same as in the unpolarized case |L2 
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where Am™ = — m" , 7711,2 being the mass of the heavy quark Qi,2 and AP qg (x) 
(2x — l)/2 the polarized g — > g splitting function. 



The ^-integrated results in [] 1] are easily recovered if one notices that 1/(1 — Zh) 

L = ]n(l-Zh,min)/0--Zh,max), 1/ *h " > L = In 4,mai/Zft,min = <~> ™ 2 ], 1/(1-^) 2 



W 2 /m 2 and l/z^ — > t>t> W 2 /m\ with the ^-integration boundaries [14 



,Q 2 ) = -[l + /\m 2 /W 2 ±vv] 



;n) 



where u 2 (x, Q 2 ) = 1 - (mj + ma) 2 /^ 2 , v 2 (x, Q 2 ) = 1 - (m x - m 2 ) 2 /W 2 and W 2 (x, Q 2 ) = 

Q 2 

— (1 — x). In the following we are only interested in neutral current charm production, 

x 

i.e., mi — m 2 — m c, v = 1 and L = L = ln(l +v)/(l — v). In this case the ^-integrated 



partonic structure functions read: 
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2.2 Electron-quark scattering 



For completeness we list the LO neutral current results for electron quark scattering [|TT|][|: 
9 BB '(x,Q 2 ) = \ J2 S^'[Aq(x,Q 2 )+Aq(x,Q 2 )] 

q=u,d,s 

gi B '(x,Q 2 ) = - <+'[Ag(x,Q 2 )-Ag(x,Q 2 )] (13) 

q=u,d,s 

gf B '(x,Q 2 ) = 2xgi B '(x,Q 2 ). 



^ote that = -R q /2 with R q defined in (Tl| 
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3 Radiative corrections in LLA 



As already mentioned, in 0(a)-LLA the structure of the QED-corrections to polarized 

n _|_ x 2 1 

DIS is the same as in the unpolarized case, because AP ee (x) = P e e(x) = • We 

L l — x j , 

therefore restrict ourselves to give the formulae that we have used in this work. For details 



we refer the reader to the literature |7|, [13 

We are using leptonic variables 

qi=p e - p e , Qi = ~q t , Vi = , xi = = — , (14 

V ■ Pe 2p-qt Syi 

which are reconstructed by the 4-momentum p' e of the observed electron and the 4- 

momenta p e , p of the initial state electron and proton. As usual S = (p + p e ) 2 is the 

center-of-mass (CMS) energy. 

To describe semi-inclusive processes we further define the Lorentz-invariant variable 

z D = z D i = (15) 
P-Qi 

which is built from the 4-momentum p^ of the produced D-meson. In the rest frame of 
the proton Zp is the energy fraction carried by the D-meson with respect to the energy 
of the boson. 

The rescaled variables for the partonic subprocess like ,e.g., 
e(p e = zxp e ) + g(p = z 2 p) — ► e(p' e = p'Jz z ) + c{p c = P D /z A ) + X read: 

v a 2 ^2 z i n 2 « M ziz 3 -l+yi 

q = Pe-Pe, Q =-q = —Qi, y = — = , 

Z 3 PPe ZiZ Z 

X = 7-—: = -j- = —, — r, S = ip e +p) 2 = Z 1 Z 2 S, (16) 

2p-g Sy z 2 {z 1 z 3 - I + yi) 

. z D yiz 3 z D 

z c = P- Pc/P ■ q = —, ■ t- = — r 

z A {ziZ 3 + yi - I) z A 

with r defined by 



yi z 3 



z\z 3 +yi-l 

The index I to denote leptonic variables will be suppressed from now on. 



(17) 
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3.1 Inclusive case 



The corrections due to initial state radiation (ISR) and final state radiation (FSR) can 
be written as M 



d 2 Aa L 
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with a the electromagnetic coupling, L e = \nQ 2 /m 2 e (m e = electron mass) and 
H(z) 
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Actq depends on the hard scattering process. 



1. Boson gluon fusion (heavy flavour production): 
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((.9i( x ,Q 2 ) from eq. ([12])) and the Jacobian 
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The integration bounds are given by: 



y 



+ 



4m 2 



1 — xy ' S^l — xy) 
l-y + xy 



4m 2 Z3 
H -— 

q 2 *i 



(24) 
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2. Light flavour excitation: 

In the case of the scattering off the light (massless) flavours (u, d, s) the partonic cross 
section is proportional to 5(1 — x) and one easily finds (of course Ag — > Aq, q — u,d, s) 

4na 2 y 
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Sx 2 y 2 yz\ 



zixy 



X N ( y(2 - y)z 2 G x (z 2 , Q 2 ) + fz 2 G 3 (z 2 , Q 2 ) + (1 - y)G 4 (z 2 , Q 2 

(25) 



- and the structure functions G{ from eq. (|8|) and (0). The inte- 
gration bounds can be obtained from the above z™ n by setting m c — > 0. 



3.2 ^-differential case 



In this subsection we consider z^-differential polarized heavy flavour production. The 
cross sections for ISR and FSR have the same structure as in eq. ([Of ) and (|T9|) with a 
^-dependent function Aa 0: 



d 3 Acr* 
dxdydz D 



d 3 Agf 
dxdydz D 
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(A<t (1, z 3 ; z D ) - Aa (l, 1; z D )) 



d(x, y, z c ) d 3 Aa 
o(x, y, zc) dxdydz c 
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The partonic cross section is given by 
d 3 Aa 4ira 2 



dxdydz c Sx 2 y 2 



■X N y(2 - y)xGi(x, Q 2 , z c ) + y 2 xG 3 (x, Q 2 , z c 
+ (l-y)G 4 (x,Q 2 ,z c 
inserting the ^-differential structure functions from eq. ( |T0D in 

G 1 (x,Q 2 ,z c ) = A^Cf B ' X B(g 2 )XB'(g 2 )^'(x,Q 2 ,4) 

B,B' 

G 3 4x,Q 2 ,z c ) = Y, c ^ B '^{Q 2 )xB'{Q 2 )glf{&,Q 2 rz c ). 

B,B' 

The Jacobian in eq. ( f2~8|) can easily be calculated as 

d(x, y, z c ) d(x, y) dz c y 2 



d(x, y, z D ) d(x, y) dz D z 2 z 4 (z 1 z 3 -l + y) 2 ' 
The additional integration bounds read: 
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with r from eq. (]17j) an d z, 



zmin.max _ rmn,max 



h=c 



(x, Q 2 ) from eq. (|TT|) . 



For D c (z), we use a Peterson et al. fragmentation function (15 
D c (z) =n{z[1-z- 1 -e c /{1-z 
normalized to dzD c (z) = 1, i. e. , 



(4 - e c )V4ic - e 2 



V4:S C - s 2 



V4e c -e c 2 ) 2 



(29) 
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(e 2 -6e c + 4) \ e c 2-^1 1 1 

JM = - : — - < arctan — - + arctan — > H — in s c + 



4 - e c 
(34) 



4 Numerical results 

The production of heavy flavours will be measured in several future experiments, e.g., 
the charm upgrade of the HERMES experiment ||, the COMPASS experiment || and 



10 



a possible polarized HERA collider mode M. In the following, we discuss the radiative 
corrections to polarized heavy flavour leptoproduction in the kinematical regions of these 
experiments. 



As usual we introduce the correction factor 6 by dAa = dAa°(l + 5), i. e., S 
5 i + 5 f + . . . with 

§U( X y ) = d 2 Aa^ d 2 Aa° 
dxdy dxdy 

or, in the ^-differential case, 

sif , . d 3 Aa^ . d 3 Aa° 



dxdydzp dxdydzp 

The labels "i" and "f" stand for "initial" and "final" state radiation, dAa denotes the 
Born cross section under consideration. 

In the following figures we employ the polarized parton densities GRSV(LO) 'Stan- 
dard' [15 with a positive definite polarized gluon and GS(LO) 'C |T7] with a non-positive 



definite gluon. It should be stressed again that both of these different parametrizations 
describe present data. Furthermore, theoretical uncertainties due to different choices of 
the charm mass m c and the factorization scale fi can be sizeable not only in the cross 
sections but also in the ratio 5. In all figures the charm mass has been set to m c = 1.5 
GeV. An appropriate choice of fi is not known. In order to find a factorization scale \i 
leading to a good perturbative stability affords the knowledge of the next-to-leading or- 
der corrections to polarized heavy flavour production which have not yet been calculated. 
Thus, we use the two factorization scales /if = Am 2 c and fi 2 , = Q 2 + 4m 2 , which have been 



favoured in the unpolarized case |]T8|. For completeness, all figures have been calculated 
with A = Ajv = 1- 

4.1 HERMESC 

The detector upgrade of the HERMES experiment p| is dedicated to measure polar- 
ized leptoproduction of open charm via the scattering of longitudinally polarized elec- 
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trons (E e = 27.5 GeV) off a longitudinally polarized fixed nucleon target, i. e. , at 
y/S = 7.4 GeV. In fig. [I] [m c = 1.5 GeV, fi 2 = Q 2 + 4m 2 « 4m 2 , parton distribu- 
tions: GRSV(LO) 'Standard' (solid line), GS(LO) 'C (dashed line)] we show the radia- 
tive corrections in C(a)-LLA to this process according to eq. ([18|) (pQ|) for two values 
of Bjorken-x in dependence of y. From the threshold condition s = (p c + p 5 ) 2 > 4m 2 
follows y > 4m 2 /(S'(l — x)) =: y min ~ 0.16/(1 — x). In the range y m in < V < 0.4 the 
corrections are big and strongly depend on the parametrization of the polarized gluon. 
In the case of the GS parton distributions the correction factor 5 diverges at y ~ 0.30 for 
x = 10~ 2 and y ~ 0.35 for x = 10 _1 due to a change of sign of the Born cross section 
in the denominator. However, in this range the cross section dAa is small. In the limit 
V — > Umin dAa vanishes. In the physically relevant region y > 0.5 the corrections are 
smaller than 20% further decreasing with increasing y. Uncertainties due to the parton 
distributions are smaller than 5%. Using a factorization scale /i 2 = s leads to negligible 
changes in figure 0. 



4.2 COMPASS 



The COMPASS experiment 0] is going to measure polarized open charm muoproduction 
at an energy of y/S = 14.1 GeV (E^ = 100 GeV). The radiative corrections for muopro- 
duction can be obtained by simply replacing L e = lnQ 2 /m 2 — > lnQ 2 /m 2 = L^. Since 
nothing else is changed the corrections 5(x, y) are expected to be similar to fig. [I] but 
suppressed by the factor L^/L e < 0.5. Fig. || [m c = 1.5 GeV, fi 2 = Q 2 + Am 2 w 4m 2 , par- 
ton distribution: GRSV(LO) 'Standard'] shows the correction factor 5(Q 2 , v) (y = v/E^, 
x = Q 2 /2Mu) for the (z/, Q 2 )-range given in table 3.1 of 0. Only for v = 40 GeV and 
Q 2 > 6 GeV 2 the corrections exceed 5%. 
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4.3 Polarized HERA collider mode 



In a polarized HERA collider mode longitudinally polarized electrons are scattered off 
longitudinally polarized protons with a center-of-mass energy VS = 300 GeV. One of 
the main goals of a polarized HERA is the measurement of the spin dependent structure 
function gi(x,Q 2 ) at values of (x, Q 2 ) which are not accessible to other experiments ||. 
In the range of small x < 10~ 3 the charm contribution g1(x, Q 2 ) to g{ {x, Q 2 ) has to be 
accounted for. 

In fig. |3]the charm contribution to the radiative corrections to d 2 Aa/dxdy(ep —> eX) is 
illustrated. On the right side we show the electromagnetic (em) cross sections dAa 0,l (ep —>■ 
eX), including the u, d, s contributions from electron quark scattering (eq. (|28D ) and the 
charm contribution due to the photon gluon fusion (eq. (PO])), at \/S = 300 GeV for 
x = 10~ 2 , 10~ 3 , and 10" 4 , using the GRSV(LO) 'S' |6| (solid line) and the GS(LO) 'C' 
|17| (dashed line) parton densities and a factorization scale /z 2 = Q 2 + 4m 2 with m c =1.5 
GeV (for the charm contribution). The cross sections for FSR have similar shapes as for 
ISR. The corresponding correction factors S(x,y) can be seen on the left hand side. The 
dotted lines have been calculated without heavy (charm) quark contribution. Employing 
the GRSV parton distributions they agree with the dashed lines in the figures 11 and 
12 in 0. However, for x < 10~ 2 the corrections strongly depend on the choice of the 
parametrization of the polarized parton densities as can be seen by comparison of the solid 
with the dashed lines. This point has not been hinted at in the literature. The charm 
contribution to the radiative corrections is roughly 10% in regions where the corrections 
do not exceed 50% and even bigger otherwise. 

Similar to the unpolarized case (F£) [TP|, it would be interesting, albeit very difficult, 
to measure g\ in order to constrain the polarized gluon density Ag(x, Q 2 ) in the proton at 
small x. Fig. || shows the (inclusive) radiative corrections to polarized charm production 
via the photon gluon fusion. Z-exchange has been neglected, since g^ = g% = (inclusive 
case, see eq. (0)) and g^ 11 ~ g\. In contrast to unpolarized heavy flavour production 
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|7]] the theoretical uncertainties are not under control. The use of the two factorization 
scales (x\ = Am 2 and /i 2 , = Q 2 + Am 2 |18|] and the parton distributions GRSV(LO) 'S' and 
GS(LO) 'C' in each case results in rather different corrections, especially for x = 10~ 3 . 
This behaviour is not only due to the unknown NLO corrections but mainly due to the 
almost entirely unconstrained Ag(x,Q 2 ). 

As has been already stressed in measurements of heavy quark production processes 
(one hadron inclusive processes) are of course differential in the momentum of the (ob- 
served) heavy quark, recommending to perform the radiative corrections on the same 
differential level. Thus, we have repeated the semi-inclusive ^-differential analysis of 
our study of unpolarized heavy quark production. Contrary to the inclusive case, 

and #4 do not vanish (see eq. flltiD ) and should not be neglected. Allthough g±(x,Q 2 ) ~ 

Q 2 

Xz{Q 2 )gJ Z {x,Q 2 ) is suppressed at small Q 2 <C Mf by the factor XziQ 2 ) ~ O-^T^ TTf 

its weight relativ to g±(x, Q 2 ) in the cross section is (up to terms of the order 0{y 2 )) 

(i - y)94(x,Q 2 ) _ xz gj z ^ 



2xygl(x,Q 2 ) 2xy gl 

due to the structure of the hadronic tensor. The same holds for the radiative correc- 
tions where the corresponding factor in Aa (zi, z 3 ) reads v)9i{z2-,Q ) ^ factor 

2xygi{z2,Q 2 ) 

Xz/2xy w 0.7Sxy/2xy -1/(Q 2 + M^) is varying between approximately 0.7 (Q 2 = S = 10 5 
GeV 2 ) and 7 (Q 2 M§ (small x)). For this reason it is (a priori) not allowed to disregard 
g>4 (,i.e. , contributions from Z-exchange)0. However, in the case of ^-differential heavy 
flavour production the differences between em and ew radiative corrections turn out to be 
negligible. This could be expected because the electroweak effects from g 3 4 have to vanish 
after ^-integration (see eq. (pT2|) . In fig. |5] we show the ^-differential radiative correc- 
tions according to eq. (|26l)-(|28D after additional integration over the kinematical range 
10 GeV 2 < Q 2 < 100 GeV 2 , 0.01 < y < 0.7. The left hand side shows the correction 
factor S(zd) (initial and final state radiation), employing the GRSV(LO) 'S' |16| (solid 
line) and the GS(LO) 'C [17j (dashed line) parton distributions. Again we have chosen 

2 The importance of the electroweak contributions can be seen in fig. 12 of j6|.We can confirm this 
result and are grateful to D. Bardin for discussions on this point. 
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fj 2 = Q 2 + 4m 2 . with m c = 1.5 GeV. For D c a Peterson et al. fragmentation function 
15 1 with e = 0.15 has been taken. On the right side of fig. |5| the ZD-dependence of the 



cross sections dAa 0,t is displayed. As in the inclusive case the ISR and FSR cross sections 
have similar shapes. Apart from irrelevant regions of the phase space where the Born 
cross section vanishes the corrections amount between 20% and —40%. Applying suitable 
(depending on the parton distribution) cuts, e. g. , 0.2 < z D < 0.8 it is possible to further 
reduce the corrections as can be seen from the rhs of fig. [|. 

To conclude, at small x it is (up to now) only possible to perform calculations of 
radiative corrections to charm production modulo the choice of the parametrization of 
the polarized gluon. The situation will improve as soon as the gluon distribution becomes 
more constrained (for x < 10~ 2 ) and the NLO expressions for polarized heavy flavour 
production are calculated. 



5 Summary 

We have calculated fully massive QED corrections to polarized heavy flavour leptopro- 
duction in C(«)-LLA for HERMES, COMPASS, and HERA kinematics. The radiative 
corrections to the first two experiments are smaller than 20% and 5% respectively in rel- 
evant realms of the phase space and the theoretical uncertainties due to different choices 
of the parton distributions and the factorization scale are under control. For HERA kine- 
matics the situation changes if x < 10~ 2 . In this case all results strongly depend on the 
chosen set of parton distributions. However, for any fixed parton distribution we find the 
charm contribution to the corrections to a measurement of g^{x,Q 2 ) non-negligible (for 
x < 10~ 3 ). Finally, we have discussed radiative corrections to gl(x,Q 2 ) at a polarized 
HERA in analogy to F£(x, Q 2 ) |19|, ^(J in the unpolarized case. Besides the inclusive case, 
we have derived general ^-differential expressions for the polarized boson gluon fusion 
process. These formulae have been used to study ^-differential corrections, keeping in 
mind that heavy flavour production is a semi-inclusive process recommending to calculate 
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the corrections on the same differential level. These corrections turned out to lie between 
20% and —40% apart from regions of the phase space where the Born cross section is 
small. 
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Figure 2: Radiative corrections to longitudinally polarized heavy quark muoproduction in 
C>(a)-LLA at V^S = 14.1 GeV (COMPASS). 
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Figure 3: Right side: dAa°' l (ep — > eX) (u,d,s and heavy (charm) quark contribution) at 
VS = 300 GeV for x = 10~ 2 , 10~ 3 , and 10~ 4 , using the GRSV(LO) 'S' pj (solid line) and 
the GS(LO) 'C [0](dashed line) parton densities and a factorization scale /j 2 = Q 2 + Am 2 c 
with m c = 1.5 GeV. Left side: The corresponding correction factor S(x,y). The dotted lines 
have been calculated without heavy (charm) quark contribution. 
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Figure 4: Radiative corrections to longitudinally polarized heavy quark electroproduction in 
0(a)-LLA at VS = 300 GeV (HERA). 
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Figure 5: ^-differential radiative corrections to deep inelastic charm production (m c = 
1.5 GeV, jj 2 = Q 2 + Am 2 ) , integrated over the kinematical range 10 GeV 2 < Q 2 < 
100 GeV 2 , 0.01 < y < 0.7, using the parton distributions GRSV(LO) 'S' |I§ (solid line) and 
GS(LO) C [0] (dashed line) and a Peterson et al. fragmentation function [fKJ with e = 0.15. 
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